Observation of Feshbach resonances between ultracold Na and Rb atoms 
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We have successfully prepared an optically trapped ultracold mixture of 23 Na and 87 Rb atoms and 
studied their interspecies Feshbach resonances. Using two different spin combinations, several s-wave 
and p-wave resonances are identified by observing high inelastic losses and temperature increases 
for both species near resonant magnetic field values. The two s-wave resonances observed below 
500 G between atoms in their lowest energy levels are possible candidates for Feshbach molecule 
association. Our results are well characterized by a coupled-channel model which also refines the 
ground state interaction potentials between 23 Na and 87 Rb. This work opens up the prospect for 
preparing ultracold ensembles of ground-state bosonic NaRb molecules which are chemically stable 
and can provide strong dipolar interactions. 
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There has been a tremendous effort in producing ultra- 
cold polar molecules for more than one decade 0, Q . The 
successful preparation of fcrmionic 40 K 87 Rb molecules 
near quantum degeneracy in 2008 is a landmark and also 
provides us so far the most promising scheme in this en- 
deavor As both methods, atom cooling and Feshbach 
molecule association [3, Q are well developed techniques 
in ultracold quantum gas research, once proper transition 
frequencies for the stimulated Raman adiabatic passage 
(STIRAP) are found starting from Feshbach molecules, 
this scheme can reliably produce a variety of ground-state 
molecules [3, H, • The first set of experiments with KRb 
molecules has provided a glimpse of the rich amount of 
physics which can be studied with ultracold samples of 
polar molecules [8l-fTo| . 

However, the majority of the recent theoretical propos- 
als with dipolar interactions require a quantum degener- 
ate sample [TTJ] . For KRb, the large inelastic loss induced 
by chemical reaction makes further cooling to degeneracy 
problematic. But it is believed that applying the same 
production scheme developed for KRb to other molecular 
species where the AB+AB — > A2+B2 reaction is energet- 
ically forbidden at ultralow temperatures, Bose-Einstein 
condensates (BEC) or degenerate Fermi gases of polar 
molecules could be achieved. As it has been pointed out, 
alkali diatomic molecules NaK, NaRb, NaCs, KCs and 
RbCs satisfy this requirement [12| and they ar e ap plied 
presently by different groups around the world [13l - fl7| . 

In this work we present our results on Feshbach reso- 
nances (FRs) between 23 Na and 87 Rb atoms. We choose 
to work with this combination because NaRb molecule 
has, besides its chemical stability, a rather large perma- 
nent electric dipole moment of 3.3 Debye [l8|. Another 
advantage is that several low lying electronic potentials of 
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NaRb molecule, including the ground states x £ + and 
a 3 S + , were already studied in detail [l9l - l22j and pre- 
dictions of FRs could be made from these results [20| . 
We have observed 3 s-wave resonances in two different 
spin combinations which are within 50 G of these original 
predictions. A full coupled-channel calculation identifies 
several other resonances being of p-wave nature. This 
work opens up the possibility of producing NaRb Fesh- 
bach molecules and is an important step toward our quest 
for ground-state polar molecules. 

Although Na and Rb are two of the most popular 
atoms for quantum gas research, their optically trapped 
ultracold mixture has never been produced to our knowl- 
edge. We prepare them for the first time in the same 
single-chamber setup as previously described for our Rb 
BEC production [23| . Similar to that work, we use ultra- 
violet light induced atom desorption to enhance load- 
ing of the magneto-optical trap (MOT) for both Rb and 
Na. Because of the low Na saturated vapor pressure at 
room temperature, our Na MOT contains only several 
million atoms which is almost 50 times less than Rb. 
After optical pumping both Rb and Na atoms to their 
|/ = l,m,f = —1) hypcrfine states, they are loaded into 
a magnetic quadrupole trap. Forced microwave evapo- 
ration is then performed on Rb, and Na is sympathet- 
ically cooled by Rb. A 1070 nm laser beam with 4 W 
power and a 65 /mi beam waist is applied for the optical 
dipole trap (ODT). The focus of the ODT is displaced 
by 90 /mi from the quadrupole trap center to partially 
suppress Majorana losses. In this way, we can achieve 
temperatures below 10 /zK, albeit with low evaporation 
efficiency. The atoms are then loaded into a crossed ODT 
by superimposing a second laser beam also of 65 /urn 
beam waist to the focus of the first trapping beam. The 
two ODT beams cross each other with an angle of 62° 
in the horizontal plane. A low homogeneous magnetic 
bias field is applied and both Rb and Na atoms can be 
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FIG. 1. (Color online) Feshbach resonance at 347.8 G for spin 
channel / = l,m/ = lof both atoms, observed by loss and 
temperature increase on Na species, left side, and Rb species, 
right side. 

transferred into their lowest hypcrfine state |1,1) simul- 
taneously with a radio frequency (RF) adiabatic rapid 
passage. 

To further increase the mixture's phase-space density, 
we continue the evaporative cooling in the crossed ODT. 
At 1070 nm, the relative trap depths Uj^/U^a is 3. 
When ramping down the laser powers, evaporation hap- 
pens mainly for Na because it sees a much shallower trap. 
The heavier Rb is sympathetically cooled by Na. For 
the best result, it is thus beneficial to prepare initially 
more Na atoms than Rb in the crossed ODT. In the end, 
we typically have 1.5 xlO 5 atoms for each species. The 
sample temperature is 1.5 (1.8) /xK for Na (Rb). Atom 
numbers and temperatures are measured by standard ab- 
sorption imaging method. The final trap frequencies are 
measured to be 2ttx (280, 310, 140) Hz and 2ttx (320, 355, 
160) Hz, respectively. Correspondingly, the calculated 
peak density is 10 13 /cm 3 (3xl0 12 /cm 3 ) for Rb (Na). 

Coarse searching of FRs is performed by first ramp- 
ing the homogeneous magnetic field to a desired value 
and then applying a 2 Gauss sweep during the 500 ms 
holding time. The interspecies resonances manifest them- 
selves as number losses and temperature increases in both 
species. We also confirm experimentally that atom losses 
and temperature changes do not occur for these field val- 
ues with Na or Rb only. Once a resonance is located, 
the accurate resonance center is determined by the fast 
ramping and holding method with finer steps. The typ- 
ical hold times are 50 ms in these final determinations. 
Observations were done over weeks and recordings were 
repeated to evaluate their reproducibility. We found an 
error of ±0.25 G, deriving from long-term fluctuations in 
the setup. We calibrate the magnetic field by RF and 
microwave spectroscopy with a 0.1 G resolution on the 
precisely known Zeeman structure of Na or Rb and esti- 
mate the overall uncertainty to be ±0.3 G taking repro- 
ducibility and calibration into account. 



Two s-wave resonances at 347.8 G and 478.8 G were 
observed for Na |l,l)+Rb [1,1) mixtures and Fig. Q] 
shows such recordings as an example for the loss and 
temperature variations. On resonance, more than 60% 
of the atoms are lost for both species and the sub- 
ensemble of each species also shows more than 30% tem- 
perature increase. These effects are attributed to the 
enhanced heteronuclear three-body collisions near an in- 
terspecies Feshbach resonance. We fit the atom loss pro- 
files phenomenologically to a Gaussian function for ex- 
tracting tthe resonance centers B exp and the full-width- 
half- maximum(FWHM) AB exp . The uncertainty from 
the profile fit is insignificant compared to the overall un- 
certainty. The s-wave resonances have measured widths 
of about 1 G. Another s-wave resonance is observed for 
Na |l,-l)+Rb |1, -1) atoms at 899.8 G with a reso- 
nance width of only 0.4 G. As the derived widths are sig- 
nificantly influenced by three-body collisions, theoretical 
modeling is necessary to determine the two-body cou- 
pling strength between the open and closed channels. 

Based on the previously measured NaRb X 1 E + and 
o 3 E + state potentials [2(|, a full coupled-channel model 
is constructed to fit the measured FRs. The Hamilto- 
nian is set up as for several other heteronuclear alkali 
molecules [lij], and includes the Fermi contact interaction 
for Na and Rb and the atomic Zeeman effect. Initially, 
the atomic parameters for these interactions are directly 
incorporated as compiled by Arimondo et al [25[. With 
this first approach, we already can assign the FRs with- 
out ambiguity, as given in Table HI The column "Type" 
identifies the s- or p-wave character, the projection M 
of the total angular momentum and the total angular 
moment / of the coupled close channel to which the 
resonance correlates at zero field. The column "Closed 
channel" gives the vibrational assignment and the atomic 
asymptote (f/va±fRb) from which the vibrational levels 
are counted, i.e. -1 represents the first level below the 
specified asymptote. 

The model also identifies several other observed fea- 
tures to be p- waves in nature. Detailed scans reveal their 
structures related to different projections m; of the ro- 
tational angular momentum I onto the axis defined by 
the magnetic field and the resulting manifold of M = m / 
(prepared input channel) +m; . For the resonance around 
954.5 G observed with the Na|l,-1)+Rb |1,-1) spin 
combination, a 0.3 G splitting is clearly resolved as shown 
in Fig. [2] (a) by loss of Na atoms when the B field 
is scanned near these resonances. Similar loss profiles 
for Rb atoms as well as heating of both species are also 
observed. In Fig. [2] (b) p-wavc resonances for the Na 
l,l)+Rb |1,1) entrance channel are shown. The left 
two features are centered at 284.1 and 284.2 G and are 
barely resolvable due to the limited magnetic field reso- 
lution and possibly to the 1.8 fiK sample temperature. 
In Table U the type of resonance is specified by m/ of the 
closed channel, resulting from the relation M = writ + mi . 
The manifold of possible mf is determined by the corre- 
sponding / at zero field. 
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TABLE I. Observed 23 Na + 87 Rb interspecies FRs below 1000 G for two different spin combinations with the assignments 
from coupled-channel calculations. The column "I obs-cal" gives the differences between observed and calculated resonance 
positions assuming two-body collisions applying fixed hyperfine interaction. The column "II obs-cal" gives the differences for 
independent hyperfine interaction as modeled by eq. ©. The columns "Type" and "Closed channel" are explained in the text. 
Theoretical near resonance scattering length profiles are calculated for s-waves at 1 nK and fitted to eq. ((3| by the parameters 
B reB , A and a& 9 . 
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I obs-cal 


II obs-cal 
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Closed channel 
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ab 9 (a ) 


Na|l, 1 > + Rb|l, 1 > 


284.1(3) 


0.12(4) 


-0.096 


-0.051 


p M=l and 2 m f =2 


v=-2(2+l) 








m/=2 


284.2(3) 


0.14(2) 


-0.007 


-0.024 


p M=3 m f =2 












284.9(3) 


0.04(1) 


-0.815 


0.112 


p M=2 m f = l 












285.1(3) 


0.06(1) 


-0.884 


0.068 


p M=l m f = l 












347.8(3) 


0.95(17) 


-0.239 


0.009 


s M=2 / = 3 




347.75 


-4.89 


66.77 




478.8(3) 


0.99(7) 


0.066 


-0.039 


s M=2 / = 2 




478.79 


-3.80 


66.77 


Na|l, -1 > + Rb|l, -1 > 


899.8(3) 


0.45(5) 


-0.173 


0.022 


s M=-2 / = 2 


v=-l(2+l) 


899.82 


0.333 


75.91 


m/=-2 


954.2(3) 


0.12(1) 


-0.342 


-0.059 


p M=-2 m f =-2 
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-0.026 


p M=-3 and -1 m f =-2 












954.0 



954.5 955.0 



S 1.0 



0.5 




284.0 284.5 285.0 

Magnetic Field (G) 



FIG. 2. (Color online) p-wave FRs for two different entrance 
channels, (a) Na and Rb atoms are both in the |1, —1) hyper- 
fine Zeeman level, (b) Na and Rb atoms are both in the |1, 1) 
hyperfine Zeeman level. Solid circles are measured remaining 
Na atoms at different applied magnetic fields. Few error bars 
shown reflect statistical spread of several shots at the same 
field. Solid lines are fits to multiple-peaks Gaussian functions 
for extracting the resonance centers. 



For describing the p-wave splitting we have to include 
in the Hamiltonian the effective spin-spin interaction (26j 
resulting from the pure dipole-dipole coupling of the two 
electron spins and from second order spin-orbit interac- 
tion. We use a simple functional form for the effective 
molecular parameter A(i?) describing the magnitude of 
this interaction as done e.g. in [26j : 



A 



1 



a S o exp (-W?) 



(1) 



Here a is the fine structure constant, aso gives the mag- 
nitude of the spin-orbit contribution and b is set for the 



slope of the exponential to 0.7196 a^ 1 which equals to 
the choice for Rb in [27| assuming that in NaRb the spin- 
orbit contribution will be determined by Rb. We fit the 
parameter aso, which yields ago = -1-0(4) a^ 3 from the 
observed splittings. The contribution by second order 
spin-orbit coupling is small but significant at long range. 

Precise determination of FRs relies on accurate po- 
tential curves for the and a 3 S + states of 
NaRb. Therefore, both potentials are constructed in a 
power series of the internuclear separation R, with the 
parametrization obtained from earlier spectroscopy work 
[2fj |. In a series of iterations between coupled-channel 
calculations for the FRs and single channel calculations 
of rovibrational levels of the X 1 S + and a 3 S + states, the 
potentials are improved by fitting the underlying coeffi- 
cients [28| to match the locations of the FRs and simul- 
taneously about 10000 rovibrational transitions from the 
Fourier-transform molecular spectroscopy [20| . The cal- 
culated singlet and triplet levels responsible for the FRs 
are added as observables to the set of transitions with 
proper weights for joining the two fits. 

In Table U we summarize all observed interspecies 
Na+Rb FRs together with their assignments from the 
coupled-channel calculations. The deviations between 
experiment and theory, shown in column "I obs-cal" , are 
mostly within the experimental uncertainties. But the 
group of p-wave resonances around 285 G shows pecu- 
liarities within its five members (two resonances are al- 
most degenerate, as given in the table). Looking into 
the detailed level scheme we find that this structure is 
a result of a strong mixing of the level / = 2, I = 1 
and v' = —1 of the asymptote Na |1) + Rb |1) with 
/ = 3, I = 1 and v' = —2 of the asymptote Na |2) + Rb 
|1). These levels are almost degenerate for to/ = 2 at 
a field of 285 G. Thus the structure is very sensitive to 
the energy contributions by the potentials, the effective 
spin-spin coupling and the actual hyperfine interaction. 
Assuming that the potentials and the effective spin-spin 
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coupling arc already quit well fixed, we studied the in- 
fluence of a i?-dependent hyperfine interaction. Such an 
effort is justified, because the two almost degenerate vi- 
brational levels test such possible i?-function differently. 
Such an effect has been extracted in other cases like Na 
[H] and Rb through observable shifts of the FRs or 
bound level energies. 

The initially constant atomic parameter a A of the 
Fermi-contact interaction of atom A is set up as a func- 
tion of R: 

a A (R) = a AM (l + c f e- R ' dR + d / 1 ~ ^P^' ^j , 

(2) 

which represents a changeover function with width dR to 
go from the atomic value aA,hf at large R to a molecular 
value <XM,hf(l + Cf) at small R. The additional expo- 
nent (n + 1) was chosen to avoid the divergence by the 
last term for R — > 0. We use this changeover function 
for both atoms, Na and Rb. This type of function was 
found by several trials on the fit and dR was chosen fairly 
large dR = 2.0ao and the exponent set to n = 3. The 
parameter values c/ and df are determined in the fit to 
be c/=-55.2 and df= 27.3-Oq. These values give correc- 
tions on the hyperfine coupling in the order of few % at 
long range. At short range the extrapolated magnitude 
becomes too large, but it is averaged out by the fast oscil- 
lations of the vibrational wavefunction for states close to 
the atom pair asymptote. Thus this approach only mod- 
els the hyperfine interaction in the energy regime needed 
for describing the accidental degeneracy mentioned above 
and should not be applied to deeply bound states. 

The new fit is given by deviations in column "II obs- 
cal" of Tabic Q] and shows values smaller than the exper- 
imental uncertainty. We see that not only the modeling 
for the p-wave structures is improved, but also the s- 
wave resonances are better described. We should note 
that the effective number of variables for the FRs is now 
6 compared to only 9 observations including the detailed 
structure of the p-wave resonances. This fact bears the 
risk of over fitting the observations. But the 3 s-wave 
resonances in connection with the positions of the p-wave 
groups determine already the potentials and the effective 
spin-spin contribution, comparing the results in column 
"I obs-cal" of the table. We believe that the overall im- 
provement in the fit is a clear indication that the hyper- 
fine interaction needs the extension in some form, e.g. 
the chosen one in cq. |(5J|, but which is certainly not a 



unique choice. More experimental data are desirable to 
strengthen this point, e.g. precise measurements up to 
few MHz on deeply bound states like in the case of Rb 2 

In the table we also give the resonance position B res , 
the width A and the background scattering length ab g 
for characterizing the scattering length as a function of 
the B field: 



using as many terms as s-wave resonances are present 
in this entrance channel. B res (uncertainty ±0.3G) of 
these profiles are not exactly the calculated ones in the 
fit because they were derived for a 1 nK temperature 
and the fit uses the peak of the elastic two-body rate. 
The derived a^ g listed in Table U are valid in the vicinity 
of the resonances and should be accurate to 2%. The 
scattering lengths of the uncoupled states become 106.74 
do for X 1 ^ and 68.62 ao for a 3 E + with similar accuracy. 

In conclusion, we have observed and analyzed FRs be- 
tween 23 Na and 87 Rb atoms. Two of the identified s- 
wave resonances with both atoms in their lowest energy 
levels provide favorable experimental conditions for pro- 
ducing Feshbach molecules with different closed channel 
characters (triplet (singlet) for the 347.8 G (478.8 G) res- 
onance) . These weakly bound molecules could be used in 
the future in our quest for ground-state NaRb molecules. 
Having Feshbach molecules with different triplet and sin- 
glet characters should enhance our chance of finding an 
excited-state level for an efficient STIRAP using either 
a pure singlet or a singlet /triplet mixed excited state to 
reach the lowest rovibrational level in the .X^E+state. 

For accelerating future studies on NaRb we calculated 
the expected s-wave FRs of other channels for f^ a = 1 
and /J£ = 1 up to 1000 G: 

1 1, 0> + |1,1) : 484 G and 717 G 

|1,1) + |1,0) : 388 G and 519 G (i) and 796 G (i) 

1 1, 0> + |1,0) : 411 G and 526 G and 763 G (i) 

|1,-1) + |1,1) : 517 G and 727 G 

jl, 1} + |1, -1) : 576 G (i) and 906 G (i) 
The additional label (i) indicates resonances with sig- 
nificant inelastic two-body losses. The accuracy of these 
predictions is limited by the digital step size in the calcu- 
lation to 1 G. Similar calculations can be made for other 
channels and even for Na 8o Rb knowing the full potentials 
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